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Abstract: A gas cloud near a supernova remnant (SNR) provides a target for pp-collisions leading 
to subsequent 7-ray emission through vr'^-decay. The assumption of a power-law ambient spectrum of 
accelerated particles with index near -2 is usually built into models predicting the spectra of very-high 
energy (VHE) 7-ray emission from SNRs. However, if the gas cloud is located at some distance from 
the SNR shock, this assumption is not necessarily correct. In this case, the particles which interact with 
the cloud are those leaking from the shock and their spectrum is approximately monoenergetic with the 
injection energy gradually decreasing as the SNR ages. In the GLAST energy range the 7-ray spectrum 
resulting from particle interactions with the gas cloud will be flatter than expected, with the cutoff defined 
by the pion momentum distribution in the laboratory frame. We evaluate the flux of particles escaping 
from a SNR shock and apply the results to the VHE diffuse emission detected by the HESS at the Galactic 
centre. 



Introduction 

SNRs are believed to be the primary sources of 
cosmic rays (CR) in the Galaxy. Observations of 
X-ray [12] and 7-ray emission [3, 4] from SNR 
shocks reveal the presence of energetic particles, 
thus testifying to efficient acceleration processes. 
Acceleration of particles in collisionless shocks 
is a matter of intensive research in conjunction 
with the problem of CR origin [7, 6, 11]. Cur- 
rent models include nonlinear effects (e.g., [5]) 
and treat particle acceleration using hydrodynamic 
codes (e.g., [9]). The predicted spectrum of accel- 
erated particles has a power-law form in rigidity 
with index which may slightly vary around -2. 

The VHE 7-ray emission from shell-type SNRs 
has been modelled using leptonic (inverse Comp- 
ton - IC) and hadronic (7r°-decay) scenarios. The 
leptonic scenario fits the broad-band spectrum of 
a SNR assuming a pool of accelerated electrons IC 
scattering off the interstellar radiation field produc- 
ing VHE 7-rays while the magnetic field and elec- 
tron spectrum cut-off are tuned to fit the radio and 
X-ray data (e.g., [13, 18]). The hadronic model 



fits the VHE 7-ray spectrum assuming a beam of 
accelerated protons hits a target, such as a nearby 
molecular cloud [1, 14]. The latter, if definitively 
proven, would be the first experimental evidence of 
proton acceleration in SNRs. 

The assumption of a power-law ambient spectrum 
of accelerated particles with index near -2 is usu- 
ally built into the models predicting the spectra of 
VHE 7-ray emission from SNRs. However, if a 
molecular cloud is located at some distance from 
the shock, the particles which interact with the 
cloud are those leaking from the shock and their 
spectrum is different [10]. In a toy model, the 
shock accelerates particles until the highest pos- 
sible energy is reached. At this point the shock 
cannot confine the particles any longer and they 
escape into the interstellar medium (ISM). The en- 
ergy spectrum of these particles will be monoener- 
getic with the injection energy gradually decreas- 
ing as the SNR ages. The 7-ray spectrum result- 
ing from the particles interacting with a gas cloud 
will be essentially flatter than expected. The dif- 
fuse emission detected at the Galactic centre by the 
HESS [2] may be of this sort. 
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Particle acceleration in SNR shock 



We use the steady-state diffusion-convection equa- 
tion 

3dx dp' 



where f{x,p) is the isotropic (in the local fluid 
frame) part of the particle distribution, and p is 
the particle momentum in mc units. We use a 
planar geometry and assume that the gaseous dis- 
continuity (sub-shock) is located at a; = while 
the shock propagates in the positive x-direction. 
The flow velocity in the shock frame can be repre- 
sented as V{x) = —u{x) where the (positive) flow 
speed u{x) jumps from U2 = w(0— ) downstream 
to 7io = u{0+) > U2 across the sub- shock and then 
gradually increases up to ui = u{+oo) > uq. The 
particle density is assumed to vanish far upstream 
(/ 0, X ^ oo), while the only bounded solution 
downstream is f{x,p) = fo{p) = /(0,p). We as- 
sume Bohm diffusion k{p) = Kp^ /^l +p^. The 
constant K depends on the SB / B level of the mag- 
netohydrodynamic turbulence that scatters parti- 
cles in pitch angle. However, it can be rescaled out 
of eq. (1) since we are not interested in the shock 
structure here. 

To include the back-reaction of accelerated parti- 
cles on the plasma flow the following equations 
are used in a quasi- stationary acceleration regime: 
(i) the conservation of the momentum flux in the 
smooth part of the shock transition (CR-precursor, 
x>Q) Pc+pu^ = piu\, where Pc is the CR pres- 
sure (CRs escape after crossing p = pmax), (ii) the 
continuity equation pu — piui, (iii) the Rankine- 
Hugoniot relations for the sub-shock strength = 
uo/u2 = (7 -t- 1) / (7 - 1 2i?T+iM-2) where 
M is the Mach number at a; = cxd, the precursor 
compression R = u\/uq, and 7 is the adiabatic 
index of the plasma. 

These equations self-consistently describe the par- 
ticle spectrum and flow structure. An efficient so- 
lution method is to reduce this system to one in- 
tegral equation [16]. A key dependent variable is 
an integral transform of the flow profile u{x?) with a 
kernel suggested by an exact (in the limit M ^ 00, 
Pmax — > 00) asymptotic solution of the system 
which has the foUowing form 



f{x,p) = /o(p)exp 



3k 



(2) 



where the flow potential ^ = u{x')dx' and the 
spectral index at the sub-shock and downstream 
q{p) — — dln/o/dlnp. The integral transform gen- 
erates the "spectral function" of the flow velocity 
U {p) instead of the flow velocity u {x) as foUows 



1 f° 
U{p) = - / 



exp 



g(p) 
3k(p) 



du{^), (3) 



and is related to q{p) as q{p) = dlnU/dlnp + 
S/[rsRU{p)] + 3. The equation for U {p) is 
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where t — \np, to,i = Inpo.i- The injection param- 
eter p = (47r/3)(TOC^/piuf )pq/o(po) is related to 
R: 



v = Kpo{l- R' 



Rr, 



dt' 



S Jto 



U{t'] 



-1 
(5) 



Our main goal is to calculate the flux of escaping 
particles. Since they leave the accelerator through 
the p = Pmax boundary and the solution normali- 
sation is set by the particle distribution atp = pi^j 
(injection rate), the flux will depend on the entire 
solution between pi^j and Pmax- For typical SNR 
conditions, this interval spans seven orders of mag- 
nitude or more and even small errors in the slope 
along the spectrum will result in a failure of the 
flux determination. To avoid this, we first compare 
our analytic solution with a numerical one [8]. The 
result of the comparison is illustrated in Fig. 1. To 
compare we substitute the injection rate ( i.e., the 
height of the spectrum atp = pinj) indicated in the 
Monte Carlo (MC) simulations. To calculate the 
escaping flux we use the time dependent version of 
eq. (1) using g =p^f in place of / 

dg d f , ^^g\ Idu dg 

-it^d-x\:'^<'^ix) = -,Tx%' 

and calculate the temporal variation of the total 
number of particles in the region .t > assum- 
ing Pmax — const (scc [15] for a discussion of 
this assumption, along with our assumption about 
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Figure 1: Comparison of our analytic solution with 
MC simulation [8]. Below Pi„j, it is a thermal 
Maxwellian; MC simulation parameters: M=128, 
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Figure 2: Bifurcation diagram showing injection 
rate v (which is actually a control parameter) ver- 
sus r. The insert shows the region of the "phase 
transition." The upper curve shows the particle es- 
caping flux, normalised to the plasma flux far up- 
stream and multiphed by the factor 10^, for clarity. 



monoenergetic escape and for further references). 
Note that Pmax(t) is adopted in [19, 10] 
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dx 



Therefore, from eq. (6) we obtain 



dt 



— 47rl62 



•'Pinj 



dp 



P 



(7) 



In this particle balance equation the term contain- 
ing .g(pinj) is simply due to particle injection at 
P = Pinj, while the last term is the particle convec- 
tive flux downstream. The remaining term is the 
flux of escaping particles 



- 3 
4n 



du {x) g (pmax, X) 



max max ), (8) 

where we have used definition of the spectral func- 
tion U {p) eq. (3), and the solution, eq. (2). Using 
the plasma number density rii = pi/nip one can 
obtain the normalised flux as 

Qn ^ ^uf pfnax/o (Pmax) U (Pmax) 
mUl C2 pf„j/o(Pinj) 

This quantity is plotted in Fig. 2 as a function of 
compression r. Note that the injection rate ob- 
tained from MC simulations corresponds to nearly 



the maximum escaping flux possible, where the 
function r {u) saturates. This saturation is related 
to the sub-shock smearing, i?'''+^^Af^. We be- 
lieve that the MC method overestimates the in- 
jection due to the lack of the feedback from par- 
ticle driven turbulence. It was argued [16] that 
self-regulation of acceleration, which includes but 
is not limited to the particle trapping by self- 
generated waves, leads to a significant reduction of 
the injection rate which is maintained at nearly the 
critical level where the shock compression r (i^) 
rises sharply. Fig. 2. For the calculations shown in 
Fig. 2, this approach indeed reduces the injection 
by an order of magnitude compared to MC results. 
Detailed calculations of the escaping flux will be 
presented elsewhere. Here we note that the nor- 
maUsed flux Qn/n\U\ does not depend strongly 
on M saturating (when M;>100) at Qn/niUi w 
10~^/Pi„ax for ui=0.005c. This scahng is the re- 
sult of energy requirement, since Qe = cpms^^Qn 
and Qe is the escaping energy flux which is con- 
strained by the available mechanical energy flux 
and the equipartition condition. The latter means 
that for the injection rate close to the critical, about 
half of the shock energy goes into the CR. 

Calculations 

Diffuse emission has been detected at the Galac- 
tic centre by the HESS [2]. Its intensity correlates 
with the gas density as traced by the CS emission 
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Figure 3: The spectrum of 7-rays from outer (mo- 
noenergetic protons of 25 TeV, blue-dots) and inner 
clouds (power-law with index —2.29, red-dashes); 
normalisations are arbitrary. Solid line shows the 
total spectrum. Data: HESS [2]. 



up to the >200 pc from the central source, while 
the SNR shell is about half this size, ^100 pc, as- 
suming a SNR age tsnr'^IO kyr and shock speed 
10'* km s^^. Clearly, some part of this emission 
^100 pc may be produced by particles which left 
the SNR shock r< 1 .5 kyr, if the diffusion coeffi- 
cient in the ISM Dxx{lO TeV)^l kpc^ Myr"\ 
i.e. relatively recently r ^ tsnr - 
The HESS data are obtained for the combined 
emission from all gas clouds in the region. How- 
ever, the model predicts different VHE 7-ray spec- 
tra from clouds at different distances from the 
SNR. The most distant clouds should exhibit a flat- 
ter spectrum with index close to -1; closer to the 
SNR shell the spectrum may steepen and become a 
regular power-law with index about -2 for clouds 
located at or inside the shell. To evaluate the feasi- 
bility of such scenario, we calculate the spectrum 
of VHE 7-rays assuming two components: one for 
the clouds outside of the shell (a monoenergetic 
ambient spectrum of protons) and another for the 
clouds inside the shell with a power-law index sim- 
ilar to the central source -2.29. The spectrum of 
7-rays can be evaluated using the scaling approxi- 
mation [20] as described in [17]. 

To illustrate the idea. Fig. 3 shows the calculated 
two-component spectrum from the Galactic centre 
together with the HESS data. Reasonable agree- 
ment with the data can be obtained for a proton in- 
jection energy '--'25 TeV, consistent with the SNR 



age of '--^10 kyr In the GLAST energy range the 
spectrum from the outer clouds differs significantly 
from what is expected for the usual hadronic sce- 
nario: it has slope close to -1. Gabici and Aharo- 
nian [10] came to a similar conclusion. 
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